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Abstract 

The dynamic behavior of a wave rotor topped 
turboshaft engine is examined using a numerical 
simulation. The simulation utilizes an explicit, one- 
dimensional, multi-passage, CFD based wave rotor code 
in combination with an implicit, one-dimensional, 
component level dynamic engine simulation code. 
Transient responses to rapid fuel flow rate changes and 
compressor inlet pressure changes are simulated and 
compared with those of a similarly sized, untopped, 
turboshaft engine. Results indicate that the wave rotor 
topped engine responds in a stable, and rapid manner. 
Furthermore, during certain transient operations, the 
wave rotor actually tends to enhance engine stability. In 
particular, there is no tendency toward surge in the 
compressor of the wave rotor topped engine during 
rapid acceleration. In fact, the compressor actually 
moves slightly away from the surge line during this 
transient. This behavior is precisely the opposite to that 
of an untopped engine. The simulation is described. 
Issues associated with integrating CFD and component 
level codes are discussed. Results from several 
transient simulations are presented and discussed. 

Introduction 

Wave rotors represent a promising technology for use 
as high pressure, high temperature topping cycles in gas 
turbine engines. 1 ’ 3 Predicted steady state performance 
benefits of wave rotors topping small, current 
technology turboshaft engines indicate a reduction in 
specific fuel consumption of up to 16% from the 
untopped engine, with a similar increase in engine 
power. 2 Predicted rotor metal temperatures are well 
below the combustor discharge temperature, and the 
rotor rotational speed is considerably less than that of a 
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conventional turbomachinery spool performing the 
same topping function. 

Numerical simulations of wave rotors indicate that off- 
design performance is good and that transient response 
to disturbances is stable and rapid compared with 
conventional turbomachinery. 4 Until this investigation 
however, transient wave rotor simulations did not model 
the interactions with surrounding turbomachinery (e.g. 
upstream compressor and downstream turbines) that 
would actually occur in a topping cycle configuration. 
Full, steady state, wave rotor topped engine cycle deck 
investigations have been performed but these can only 
state whether the engine will operate at one point or 
another, they cannot answer questions regarding 
transient response, stability, or controllability. 4,5 For 
example, in conventional compressor-combustor- 
turbine engines rapid fuel flow rate changes are 
restricted due to compressor surge limit considerations. 
Is there a similar limitation for wave rotor topped 
engines? How does the system respond to compressor 
inlet perturbations? 

This paper presents the results from a numerical 
simulation developed to address these and other 
questions concerning the dynamic integration of wave 
rotors with surrounding turbomachinery. The 
simulation was created by combining a dynamic, 
component level, implicit engine code, GETRAN, 67 
with a dynamic, explicit, CFD based, wave rotor code 4 
into a single integrated software simulation. The two 
codes and their integration into a single code are briefly 
described. The results from several simulations are then 
presented and discussed. Comparisons are made 
between simulations of wave rotor topped and untopped 
gas turbine engines during similar transient operations. 

The paper concludes with some remarks regarding 
implications of the results, and possible future 
improvements to the simulation. 



Simulation Description 

Both the engine code and wave rotor code have been 
described extensively in the literature. 4 * 7 As such, only 
brief descriptions will be given below. 

The GETRAN code is a one -dimensional, component 
(or turbomachinery stage) level, generic gas turbine 
engine simulation. Due to its highly modular structure, 
it can simulate nearly any combination of conventional 
gas turbine components such as turbomachinery, 
combustors, diffusers, ducting, nozzles, etc. The one- 
dimensional nature of the code allows great flexibility 
in the design process and in the simulation of transient 
operations. Code modules, each simulating a different 
engine component, communicate boundary conditions 
by the use of coupling plena (volumes). The plena 
transfer information about total pressure, total 
temperature, mass flow, and other physical conditions 
from one engine component to the next. The plenum 
module imposes the time dependent, nonlinear 
conservation equations of mass momentum and energy 
between components. Time integration is done 
implicitly, ensuring stability for the typically stiff 
system of equations. 

The second code, heretofore called WR1D, is an 
explicit one-dimensional CFD based wave rotor 
simulation which tracks the gasdynamic processes in 
each of the passages of the wave rotor as they rotate 
past the various stationary ports. Familiarity with the 
wave rotor and its operation is assumed in this paper; 
however, for reference, a schematic drawing of a wave 
rotor is shown in Fig. 1. The WR1D code, though 
relatively simple because it is one-dimensional, has 
experimentally validated sub-models for most of the 
major loss mechanisms found in wave rotors. These 
make it capable of fairly realistic simulation. Explicit 
time integration is used in this code because of the small 
time steps required to accurately capture the 
gasdynamic processes in the passages. 

Integrating the Codes 

Because of the general component-plenum structure of 
GETRAN, it was relatively simple to incorporate the 
explicit WR1D simulation code into the implicit 
GETRAN code. GETRAN does not iterate on discrete 
flowfield points or volumes as do conventional CFD 
codes. It is the plena conditions that are iterated upon 
during a given time step until the governing (difference) 
equations are satisfied. For the WR1D code, on the 
other hand, plena information are used as boundary 
conditions and port flow information is computed at 
each of its time steps. Typical GETRAN time steps are 



Figure 1 Schematic wave rotor 

several orders of magnitude larger than those in the 
WR1D code. Thus, with the WR1D simulation 
incorporated as a GETRAN component (e.g. a callable 
subroutine module), GETRAN simply supplies the 
WR1D code with it current ‘best guess’ of plena 
conditions. The WR1D code then integrates over the 
GETRAN time step and returns the port flow 
information (e.g. mass flux, enthalpy flux, etc.) to 
GETRAN. 

The “Gattling gun” nature of the wave rotor’s individual 
passages presented the only difficulty in the merging the 
two codes. The WR1D code operates on time scales 
which resolve the rapid fluctuations in flow conditions 
as the wave rotor passages open and close to the various 
ports. Typically, these fluctuations are highly damped a 
short distance down or upstream from the rotor face in 
any given duct leading to or from the wave rotor. The 
WR1D code does not model these ducts, so the 
fluctuations are not damped. The fluctuations made it 
difficult for the GETRAN iteration to converge. As 
such, it was necessary to use a simple data smoothing 
algorithm (e.g. an averaging) on the results returned by 
the WR1D subroutine. The averaging time was chosen 
so that accuracy was maintained over the time scales of 
interest in the simulations. 

Computing Requirements 

The integrated simulation code was very CPU intensive. 
Even though a very course spacing of 10 numerical 
cells per wave rotor passage was used, the small time 
step of the WR1D code (6.25xl0 6 sec.) required 
considerable CPU usage in and of itself. The implicit 
nature of the GETRAN code requires several iterations 
per GETRAN time step (8.00x10^ sec.), and therefore 
several calls are made to the WR1D routine at each 
GETRAN time step, thereby multiplying the problem. 
Typical simulation requirements were 10 hours of CPU 
time for 1 second of real time on an SGI Power 
Challenge. 
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Figure 2 Simulation schematics, with relevant component numbering for: (a) wave 
rotor topped, and (b) untopped engines 


efficiencies, and overall compressor performances, and 
Methodology shaft moments of inertia (rotor moments were 

The methodology for this investigation was relatively estimated). The first design, Case 1, was a wave rotor 

simple. Two basic engine designs were created for topped configuration. The second design, Case 2, was a 

study based upon a small turboshaft engine. Schematic baseline configuration with no wave rotor. Details of 

layouts for the two designs are shown in Fig. 2. Both the Case 1 wave rotor design has been described in the 

designs had identical design inlet mass flow rates, literature. 4,5 The wave rotor was designed to be self- 
compressor pressure ratios, and high pressure turbine driven, that is, requiring no external drive motor or 

inlet temperatures. These are listed in Table 1. Both transmission. In the present investigation however, the 

designs also contained identical design point component rotor was assumed to run at constant speed throughout 
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Table 1 Compressor and Turbine Design 
Specifications 



Case 1 

Case 2 

Design Mass Flow, kg/s 

2.168 

2.168 

Compressor Pressure Ratio 

7.77 

7.77 

HP Turbine Pressure Ratio 

2.92 

2.92 

HP Turbine Inlet Temperature, K 

1329 

1327 

HP Turbine Efficiency 

0.86 

0.86 

LP Turbine Pressure Ratio 

3.26 

2.49 

LP Turbine Inlet Temperature, K 

1056 

1052 

LP Turbine Efficiency 

0.86 

0.86 


all of the transient operations. The power required to 
maintain constant rotor speed is very small (e.g. less 
than 0.5% of the power turbine output). 


The compressor and turbines were simulated in the 
GETRAN code by use of the global compressor and 
turbine subroutines available in the code. 6 These 
subroutines model component performance (e.g. 
compressor or turbines) using mathematical curve-fits 
to the maps. The compressor performance maps and 
turbine design point efficiencies used for the curve- 
fitting were obtained from a steady-state cycle deck 
simulation. 8 

Steady State 

Both Case 1 and 2 were run until a steady state (or 
equilibrium) condition was obtained. The equilibrium 
air mass flow through the Case 1 combustor was a 
nominal 3.35 kg/s, or 54% higher than that of Case 2. 
This was consistent with specified design point values. 
Similarly, pressures and temperature in the plena 
associated with the wave rotor ports were well matched 
to design specifications. These results indicated that the 
simulation was operating correctly, at least in steady 
state. The shaft output of the Case 1 low pressure 
power turbine was approximately 0.657 MW. The non- 
topped Case 2 power output was 0.527 MW, giving the 
Case 1 configuration an expected 24.7% increase in 
power. 

It is noted that the compressor design for Cases 1 and 2 
has a small amount of bleed flow which would normally 
be used for cooling the high pressure turbine (HPT). 
Due to the addition of the wave rotor in Case 1, the 
bleed flow no longer has enough pressure to be used for 
HPT cooling. For the purposes of this investigation, it 
was simply dumped overboard. In any practical wave 
rotor topping application, this bleed flow will have to be 
extracted from the wave rotor itself; however, this is not 
expected to significantly change the dynamic results. 



Figure 3 Fuel mass flow rate schedule for Case 1 and 
2 

Transient Simulation Results and Discussion 
Two transient operations were simulated. The first 
involved rapid change in the fuel flow rate. The second 
involved a rapid drop in the compressor inlet pressure 
brought about by augmenting the simulation 
configurations with an idealized valve upstream of the 
compressor face. 

For the results to be presented, a brief word about the 
case nomenclature is in order. As previously stated, 
Case 1 and Case 2, designate the wave rotor and control 
(untopped) engines respectively. The different modes of 
transient operation are designated by a letter following 
the case number. “F S” designates the varying fuel 
schedule operation, and “V” designates the valve 
generated compressor inlet disturbance. 

Varying Fuel Flow Rate 

The fuel flow schedule used for this transient simulation 
is shown in Fig.3. The same fuel schedule was used for 
Case 1 and 2. The design point fuel flow rate for both 
Case 1 and Case 2 was 0.049 kg/s. The simulation was 
run at steady state for 0.1 sec. The fuel flow rate (m f ) 
was then rapidly ramped down to 50% of nominal over 
a period of 0.05 sec., held at that value for 3 seconds, 
then increased back to the design point value over a 
period of 0.25 sec., and held there for the duration of 
the 6.0 second simulation. 

Time traces of pressure throughout the engine, for both 
Case 1 and Case 2 are shown in Fig. 4. Combustor exit 
and high pressure turbine inlet temperatures are shown 
in Fig. 5 (they are coincident in the Case 2 
configuration). Mass flow rates at several locations are 
shown in Fig. 6. In general, all of these figures indicate 
that engine response time was comparable for both Case 
1 and 2 and in both cases the response was stable. 
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(a) 

Figure 4 Time traces of stagnation pressure at various locations of the engine for the fuel flow transient: (a) 
wave rotor topped turboshaft engine (CASE 1); (b) non-topped engine (CASE 2) 
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Figure 5 Time traces of combustor exit and high pressure turbine inlet stagnation temperature for the fuel 
flow transient: (a) wave rotor topped turboshaft engine (CASE 1); (b) non-topped engine (CASE 2) 
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(a) (b) 

Figure 6 Time traces of air mass flow at various locations for the fuel flow transient: (a) wave rotor topped 
turboshaft engine (CASE 1); (b) non -topped engine (CASE 2) 
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(a) (b) 

Figure 7 Time traces of spool angular velocities for the fuel flow transient: (a) wave rotor topped turboshaft 
engine (CASE 1); (b) non-topped engine (CASE 2) 


Spool rotational speeds are plotted in Fig. 7 as fractions 
of the design point value. The equilibrium values for 
rotational speeds were 5340 rad./s for the high pressure 
turbine-compressor spool, and 3500 rad./s for the power 
turbine spool. It is noted that the load placed on the low 
pressure power turbine essentially followed the fuel 
flow schedule except the bounds were between 40% and 
100% of the design value, and the ramp back up from 
low to high load was delayed until 1.0 sec. after the fuel 
flow rate was increased. Note that the high pressure 
spools responded identically in both Case 1 and 2; 
however, the power turbine spools were quite different. 
For the wave rotor topped case, the power turbine was 
consistently running above the design speed. The non- 
topped power turbine spool ran below the design value 
for a significant portion of the transient. The reason for 
the difference in behavior is probably due to the fact 
that at 50% of design fuel flow rate the power turbine of 
the wave rotor-topped engine can deliver 40% of design 


load at steady state. The non-topped engine can only 
deliver 28% of design point power at the 50% fuel flow 
rate. Nonetheless, the transient result implies that the 
wave rotor topped engine configuration has a better 
available shaft power response than the untopped 
engine. 

Significant differences between the wave rotor topped 
simulation and the non-topped simulation can be seen 
when the trajectory of the compressor pressure ratio (FI) 
versus corrected mass flow (me) is plotted for each 
case. This has been done in Fig. 8. As with many 
conventional turboshaft engines, the deceleration 
response (e.g. reduced fuel flow) of the non-topped 
Case 1 compressor at first follows a line of constant 
rotational speed, away from the surge line, then moves 
along a nearly straight line until a new steady state 
operating point is reached. During acceleration, the 
response is similar except the initial trajectory is toward 
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Figure 8 Compressor pressure ratio as a function of 
corrected mass flow during the fuel flow transient 

the surge line at constant rotational speed. This 
behavior means that the engine acceleration must be 
limited in order to avoid surge or rotating stall 
instabilities. The wave rotor-topped engine, on the 
other hand, traces a completely different trajectory. 
During deceleration, the compressor moves along a 
constant speed line, slightly toward the surge line, then 
it essentially parallels the surge line to the new steady 
state operating point. During acceleration, the response 
is essentially a mirror image, with initial trajectory away 
from the speed line. Overall then , Fig. 8 indicates 
significant enhancement for the wave rotor-topped 
engine in overall compressor stability during fuel flow 
transients 

Valve Generated Compressor Inlet Disturbance 
The objective of this transient was to examine the 
response of the two engine configurations to a rapid 
change in compressor inlet pressure. The inlet pressure 
change was implemented by adding an ideal valve 
upstream of the compressor face. The simulation was 
run at steady state, with the valve fully open for 0.1 
seconds, the valve opening was then instantaneously 
reduced to 20% of its fully open area and held there for 
the duration of the 3.0 second simulation. 

Pressure and mass flow rates at various engine locations 
are shown in Figs. 9 and 10. Although the response 
times of the two configurations were similar in this case, 
the amplitudes of the Case 1 were larger than Case 2. 

Rotor speed responses are shown in Fig. 1 1 . Again, it is 
seen that the changes in the wave rotor-topped engine 
were significantly larger than the untopped engine. It is 
noted that neither case actually reaches steady state. 
The load on the power turbine was maintained at the 
design value for both Casel and Case2. As such, the 


power turbine continually decelerates throughout the 
transient. 

As in the fuel flow transient, a plot was made of the 
trajectory of the compressor pressure ratio versus 
corrected mass flow. This is shown in Fig. 12. In this 
case, the wave rotor-topped engine produced a curious 
“loop” in the trajectory which was not seen in the 
untopped engine. The trajectory for the wave rotor 
topped transient moved closer to the surge line than did 
the non-topped case. However, it can be seen that little 
surge margin is required for both topped and non- 
topped transients. The explanation for this behavior is 
unclear. 

Concluding Remarks 

In conclusion several facts need to be noted regarding 
the simulation itself. The first was that the two codes 
were successfully merged despite the mixture of 
implicit and explicit solution techniques. This was 
possible and in fact, easy, due to the modular structure 
of the GETRAN code and its use of iteration upon 
solely the plena conditions. Future mergers between 
explicit 2-D and 3-D simulation routines and GETRAN 
could be possible in the future. The one disadvantage of 
the code integration was that it was very 
computationally intensive. Further research is necessary 
to improve the computational speed of the merged code. 

The wave rotor topped engine appeared robust and 
quick to respond to transient operations. In particular, 
the wave rotor topped engine demonstrated little 
tendency toward surge during rapid fuel flow changes. 

During rapid inlet pressure perturbation transients, the 
compressor of the wave rotor topped engine moved 
closer to the surge line than the non-topped compressor, 
but was well within the surge margin. 

Computed steady state performance of the wave rotor 
topped engine compared favorably with the results from 
previous studies using cycle deck calculations and wave 
rotor maps. 1 
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(a) (b) 

Figure 9 Time traces of Compressor inlet and exit stagnation pressures for the valve generated compressor 
inlet disturbance transient: (a) wave rotor topped turboshaft engine (CASE 1); (b) non-topped engine (CASE 
2 ). 
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(a) (b) 

Figure 10 Time traces of air mass flow at various locations of the engine for the valve generated compressor 
inlet disturbance transient: (a) wave rotor topped turboshaft engine (CASE 1); (b) non-topped engine (CASE 
2 ). 
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Figure 11 Time traces of spool angular velocities for the valve generated compressor inlet disturbance 
transient: (a) wave rotor topped turboshaft engine (CASE 1); (b) non-topped engine (CASE 2). 
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Figure 12 Compressor pressure ratio as a function of 
corrected mass flow during the valve generated 
compressor inlet disturbance transient. 
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